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ABSTRACT 

We present surface photometry and stellar kinematics of NGC 2950, which 
is a nearby and undisturbed SBO galaxy hosting two nested stellar bars. We use 
the Tremaine- Weinberg method to measure the pattern speed of the primary 
bar. This also permits us to establish directly and for the first time that the two 
nested bars are rotating with different pattern speeds, and in particular that the 
rotation frequency of the secondary bar is higher than that of the primary one. 



Subject headings: galaxies: individual (NGC 2950) — galaxies: kinematics and 
dynamics — galaxies: elliptical and lenticular, cD — galaxies: photometry - 
galaxies: structure 



1. Introduction 

Large-scale bars are present in some 2/3 of all disk galaxies (Knapen et al 2000; 
Eskridge et al. 2000). Secondary stellar bars within large-scale bars are also common, 
occurring in about 1/3 of barred galaxies (Laine et al. 2002; Erwin & Sparke 2002). Interest 
in secondary stellar bars is motivated by the hypothesis that they are a mechanism for 
driving gas to small radii to feed the supermassive black holes powering active galactic 
nuclei (e.g. Regan & Mulchaey 1999). However, the efficiency of such transport is uncertain 
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because of the complete lack of knowledge of VL P and Q s , the pattern speeds of the primary 
and secondary bars. Whereas a number of pattern speeds of large-scale bars have been 
measured (see Aguerri et al. 2003, hereafter ADC03 and references therein), no such 
measurements in nested systems have been performed yet. The presence of nested bars with 
different pattern speeds has been inferred largely on the basis of their apparently random 
relative orientations (Buta & Crocker 1993). The possibility that Q s > Q p is supported 
by the simulations of Rautiainen et al. (2002), which showed that secondary bars can 
naturally form and survive for more than few rotation periods in pure stellar disks. The 
morphological characteristics of these systems are suggestive of stars in the secondary bars 
oscillating about the loop orbits studied by Maciejewski & Sparke (2000) for models with 
Q s > Q p . However, simulations have also found other possibilities (Friedli & Martinet 
1993), including cases where two stellar bars counter- rotate (Sellwood & Merritt 1994; 
Friedli 1996). 

A model-independent method for measuring pattern speeds is the Tremaine- Weinberg 
method (Tremaine & Weinberg 1984, TW hereafter). This gives the pattern speed f2 of a 
single bar as 

Xnsini = V, (1) 

where X and V are luminosity-weighted average position and line-of-sight velocity measured 
parallel to the major axis of the galaxy disk, and i is the disk inclination. Long-slit spectra 
parallel to the disk major-axis can measure all the quantities needed by Eqn. 1 provided 
that the galaxy is free of extinction. If several parallel slits at different offsets Y relative to 
the major axis are available for a galaxy, fisini can be obtained as the slope of a plot of V 
versus X. In a double barred galaxy (hereafter S2B), for a slit passing through both bars, 
Eqn. 1 is modified to 

(X p Q p + X s Q s )smi = V (2) 

provided the two bars are rigidly rotating through each other, i.e. that the total surface 
density is described by S(i?, 0) = E p (i2, — Vt p t) + T> S (R, <p — Vl s t). Eqn. 2 is then a 
consequence of the linearity of the continuity equation. Eqn. 2 can be solved for fl s by first 
measuring Q p as in Eqn. 1 with slits which avoid the secondary bar, and then modeling to 
obtain X s from the observed X = X s + X p . However, the two bars are not likely to rotate 
rigidly through each other when fl p ^ Q s (Louis & Gerhard 1988; Maciejewski & Sparke 
2000; Rautiainen et al. 2002), requiring additional modeling to obtain Q s . Nonetheless, 
when Q p = Q s , Eqn. 2 reduces to Eqn. 1 and is satisfied exactly. Thus testing whether 
Q s = Q p does not require any assumptions. 
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2. NGC 2950 

An ideal target for this purpose is the S2B NGC 2950, which is a large (2'.7 x 1'.8 [de 
Vaucouleurs et al. 1991, hereafter RC3]) and bright (B? = 11.84 [RC3]) early-type barred 
galaxy. NGC 2950 is classified RSBO(r) and its total absolute magnitude is M% T = —20.03 
(RC3) adopting a distance of 23.3 Mpc (Tully 1988). The presence of a secondary stellar 
bar has been discussed extensively by Wozniak et al. (1995), Friedli et al. (1996), and 
Erwin & Sparke (2002) on the basis of ground-based and Hubble Space Telescope images in 
both optical and near-infrared bandpasses. The secondary bar of NGC 2950 is typical of 
those in the samples of Laine et al. (2002) and of Erwin & Sparke (2002). NGC 2950 meets 
all the requirements for the TW analysis: it has an intermediate inclination, both bars have 
intermediate position angles (hereafter PA) between the major and minor axes of the disk 
and no evidence of spirals, patchy dust or significant companions. 

3. Surface photometry 

The photometric observations of NGC 2950 were carried out at the 1-m Jacobus 
Kapteyn Telescope on December 27-28, 2000. We took multiple exposures in the Harris B 
(4 x 1200 s), V (3 x 480 s), and / (18 x 150 s) bandpasses using the SITe2 2048 x 2048 
CCD. This camera has a scale of 0'.'33 pixel -1 , yielding an unvignetted field of view of 
~ 10' x 10'. The seeing FWHM was w l'.'O. The data reduction has been carried out 
using standard IRAF tasks as in Debattista et al. (2002, hereafter DCA02). Images were 
bias subtracted, flatfielded, cleaned of cosmic rays, and corrected of bad pixels. The 
sky-background level was removed by fitting a second-order polynomial to the regions free of 
sources. Photometric calibration, using standard stars, included corrections for atmospheric 
and Galactic extinction, and for color as in DCA02. 

The radial profiles of surface brightness, ellipticity and PA were obtained by fitting 
elliptical isophotes with the IRAF task ELLIPSE. We first fitted ellipses allowing their centers 
to vary to test for patchy dust obscuration. We found no evidence of a varying center 
within the errors of the fits, and similar PA and ellipticity profiles for all bandpasses. Thus 
we concluded that there is little or uniform obscuration, as required for the TW method. 
The ellipse fits were then repeated with the ellipse center fixed; the resulting photometric 
profiles are plotted in Fig. 1. We interpreted the local maximum in ellipticity at r ~ 3" 
and the corresponding twist and stationary value in PA as the photometric signatures of 
the presence of a misaligned secondary bar inside the primary bar. This is confirmed by 
the analysis of B — I and V — I color maps and unsharp mask of the original frames; in 
all we find no evidence of other small-scale structures such as nuclear rings or disks, spiral 
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arms, star-forming regions, dust lanes and/or dust patches, in agreement with previous 
results (Wozniak et al. 1995; Friedli et al. 1996; Erwin & Sparke 2003). In particular, 
the structural details unveiled by the unsharp mask of the WFPC2/F814W image of the 
nucleus of NGC 2950 (Erwin & Sparke 2003) are unlike those typical of nuclear stellar disks 
(Pizzella et al. 2002). The PAs of the primary (PA P = 152?6 ± 0?4) and secondary bar 
(PA S = 91? 1 ± 0?3) were measured in the /—band image at r ~ 3" and r ~ 23", at the 
two peaks in the ellipticity profile (Fig. 1). The lengths of the primary (a p = 34'/3 ± 2'/l) 
and secondary bar (a s = 4"5 ± l'.'O) were measured using three independent methods based 
on Fourier amplitudes (Aguerri et al. 2000), Fourier and ellipse phases (DCA02), and 
a decomposition of the surface brightness profiles (Prieto et al. 2001). The inclination 
[i = 45?6 ± 1?0) and PA of the disk (PA d = 116?1 ± 1?0) were determined by averaging the 
values measured between 65" and 100" in the /—band profile (Fig. 1). 

4. Long-slit spectroscopy 

The spectroscopic observations of NGC 2950 were carried out at the 3.6-m Telescopio 
Nazionale Galileo on December 18, 2001 (run 1), March 20-22, 2002 (run 2), and March 
9-11, 2003 (run 3). The Low Resolution Spectrograph mounted the HR-V grism No. 6 
with 600 grooves mm -1 and the &!7 x 8 f 1 slit. The detector was the Loral CCD with 
2048 x 2048 pixels of 15 x 15 fxm . The wavelength range from 4660 to 6820 A was covered 
with a reciprocal dispersion of 1.054 A pixel -1 and a spatial scale of 0.275 arcsec pixel -1 . 
We obtained 4 spectra with the slit along the disk major axis (run 2 and 3), and 11 offset 
spectra with the slit parallel to it (Y = -3"1, +1"5, +2"8, ±10"1 in run 1, Y = ±5"1, ±13"4 
in run 2, and Y = — 2"5, +3"5 in run 3, Fig. 2). The exposure times were 2 x 60 min 
and 2 x 45 min for the major-axis spectra obtained in run 2 and 3, respectively, and 45 
min for all the offset spectra. Comparison lamp exposures before and/or after each object 
integration ensured accurate wavelength calibrations. Spectra of G and K giant stars served 
as kinematical templates. The seeing FWHM was ~ 1"2 in run 1, « l'/5 in run 2, and ~ 0"9 
in run 3. Using standard MIDAS routines, all the spectra were bias subtracted, flatfield 
corrected, cleaned of cosmic rays, corrected for bad pixels and wavelength calibrated as in 
DCA02. The accuracy of the wavelength rebinning (^ 2 km s -1 ) was checked by measuring 
wavelengths of the brightest night-sky emission lines. The instrumental resolution was 3.10 
A (FWHM) corresponding to a inst 80 km s _1 at 5170 A. The major-axis spectra obtained 
in the same run were co-added using the center of the stellar continuum as reference. In all 
the spectra the contribution of the sky was determined by interpolating along the outermost 
~ 20" at the edges of the slit and then subtracted. 
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5. Pattern speeds of the primary bar and secondary bar 

To measure V for each slit (Fig. 3a), we first collapsed each two-dimensional spectrum 
along its spatial direction in the wavelength range between 5060 and 5490 A, obtaining 
a one-dimensional spectrum. The value of V was then derived by fitting the resulting 
spectrum with the convolution of the spectrum of the K1III star HR 4699 and a Gaussian 
line-of-sight velocity profile by means of the Fourier Correlation Quotient method (Bender 
1990, hereafter FCQ) as done in ADC03. We estimated uncertainties by Monte Carlo 
simulations with photon, read-out and sky noise. To compute X for each slit (Fig. 3b), we 
extracted the luminosity profiles from the V — band image along the position of the slit after 
convolving the image to the seeing of the spectrum. The V— band profiles match very well 
the profiles obtained by collapsing the spectra along the wavelength direction, confirming 
that the slits were placed as intended. We used the V— band profiles to compute X because 
they are less noisy than those extracted from the spectra, particularly at large radii. We 
obtained f2 p from the values of X and V for the slits at \Y\ > 3'.'1 and at Y — 0" (the latter 
constrain only the zero point). Since the slits at \Y\ > 3'/l do not cross the secondary bar, 
we assume X = X p and V = V p for them and obtain fi p sini with a straight line fit. This 
gives tt p = 11.2 ± 2.4 km s -1 arcsec -1 (99.2 ± 21.2 km s -1 kpc -1 , Fig. 3c) The value of fi p 
does not change within errors (f2 p = 10.9 ± 2.4 km s -1 arcsec -1 ) when we exclude the slits 
at Y — — 3'/l, +3'.'5, confirming that they are not much affected by the secondary bar (Fig. 
2). 

We used the FCQ to measure the line-of-sight velocity and velocity dispersion profiles 
of the stellar component along the major axis (Fig. 4). All the major-axis spectra were 
co-added after being convolved to the same seeing. We derived the circular velocity in the 
disk region, V c = 356^9 km s -1 , after a standard correction for the asymmetric drift as in 
ADC03. Thus the corotation radius of the primary bar is D p = V c /Q p = 32"4l^ and the 
ratio Tip = D p /a p of the corotation radius to the bar semimajor axis is 1Z V = l-Ol^l (error 
intervals on D p and 1Z P are 68% confidence level and were measured with Monte Carlo 
simulations as in ADC03). We conclude that, within the errors, the primary bar of NGC 
2950 is consistent with all previous measurements of 1Z in SB0 galaxies (ADC03), which 
gives us confidence in our assumption that the signals in the outer slits are generated by 
the primary bar only. 

The photometric and kinematic integrals measured with the innermost slits (\Y\ < 2 / .'8) 
include a contribution from the secondary bar. In particular, |V| ^> |V P | for the slits at 
Y = — 2'/5, +2'.'8 if we extrapolate V p from large A straight line fit to X and V for the 
slits at |y| < 2'/ 8 has a slope (= 63.7 ±7.1 km s -1 arcsec -1 , Fig. 3c) which is different 
at better than 99% confidence level from the slope (= 8.0 ± 1.7 km s -1 arcsec -1 , Fig. 3c) 
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of the straight line fit for the primary bar. This may be because Q p ^ Q s ; however, to 
justify this conclusion, we needed to exclude the possibility of a systematic error due to 
PA errors, which X and V are quite sensitive to (Debattista 2003). We tested whether 
the difference in the slopes of the straight line fits may be due to only PA errors when 
Q p = Q s by building a model of a S2B galaxy with a single pattern speed from an N— body 
simulation of a SBO galaxy described in Debattista (2003). We rotated particles by 90°, 
rescaled their phase space coordinates by a factor of 0.2, and added them back to the 
original galaxy with various primary to secondary mass ratios to approximately match 
NGC 2950. After projecting this system as in NGC 2950, we proceeded to measure X and 
V for slits misaligned with the major axis. Even when PA errors reached ±5°, we were not 
able to produce a system which approaches the behavior of our observations. In particular, 
we were not able to produce a system in which the slopes of the integrals plotted versus \Y\ 
are larger for V and smaller for X in the region of the secondary bar than in the region 
of the main bar, as we observed in NGC 2950 (Fig. 3a,b). We therefore concluded that 
the signatures we observed could not be an artifact of any PA misalignments on two bars 
rotating with a single pattern speed. Our results, therefore, lead us to conclude, directly 
and for the first time, that Q s ^ Q p . 

However, estimating Q s is model dependent; we illustrate this by considering two 
extreme cases. First, we assumed that the secondary bar dominates at \Y\ < 2'.'8, and 
used Eqn. 1 with f2 replaced by f2 s to find f2 Sjl = 89.2 ± 9.9 km s _1 arcsec -1 (Fig. 3c). 
In the second case, we rewrote Eqn. 2 as X S (Q S — Q p ) sini = V — VL P X sini. The observed 
quantities are X and V, while fl p was measured above. To obtain X s , first we derived the 
values of X p in the region of the nuclear bar by fitting a straight line to the X values at 
\Y\ > lO'.'l (fits extending to smaller |Y| give larger \ fl s \), and obtained X s = X — X p . Then 
plotting (V — Q p X sini) versus X s , we obtain (Q s — Q p ) sini as the slope of the best fitting 
line; the result is fi Si 2 = —21.2 ±6.2 km s _1 arcsec -1 , i.e. secondary bar counter-rotating 
relative to the primary. (Note that the range from f2 Si i to f2 Sj 2 passes smoothly through 
±oo, i.e. a vertical line.) 

6. Conclusions 

We showed that the primary bar in NGC 2950 is rapidly rotating. If this is the norm 
in S2B galaxies, then it guarantees that primary bars are efficient at funnelling gas down 
to the radius of influence of secondary bars. In Fig. 4, we plot the lines of slope J\isini 
and i\ 2 sini The range of Q s is large enough that it must include the case where 1Z S ~ 1. 
However, it also includes the case where 1Z S ~ 2, which hydrodynamical simulations find 



-7- 



leads to inefficient gas transport (Maciejewski et al. 2002). 

We suggest two avenues for fruitful future work. First, since the two bars cannot be in 
exact solid body rotation (Louis & Gerhard 1988; Maciejewski & Sparke 2000; Rautiainen et 
al. 2002), a more accurate measurement of Q s will require careful modeling and comparison 
with simulations to account for such effects. Second, it may be that secondary bars oscillate 
about an orientation perpendicular to the primary bar, possibly accounting for Q s < 0. This 
can be tested by repeating our measurements on a sample of S2B galaxies. Nonetheless, 
we can confidently conclude that in NGC 2950 the two bars must have different pattern 
speeds, with the secondary bar having a larger pattern speed. 
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Fig. I. — (a) Surface brightness, //, (b) ellipticity, e, and (c) PA radial profiles of NGC 2950. 
Solid, dashed and dotted lines refer to /, V, and B— band data, respectively. The thick lines 
represent the fits to the /—band e and PA of the galaxy's disk. 
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Fig. 2. — Left panel: Large-scale image of NGC 2950 showing the primary bar and disk 
with J— band contours and slit positions overlaid. Contours are spaced at 0.5 mag arcsec -2 
and the outermost corresponds to fij = 23.0 mag arcsec -2 . The solid and dashed lines 
correspond to the position of the spectra obtained along the disk major axis and at large 
offsets > 5'.'1), respectively. For each slit position the offset, Y, is given in arcsec (we 
arbitrarily chose axes such that Y increases from the SW to the NE sides). Right panel: A 
zoom into the central region of NGC 2950 showing its secondary bar. /—band isocontours are 
spaced at 0.5 mag arcsec -2 , with the outermost corresponding to /// = 19.5 mag arcsec -2 . 
For each spectrum obtained at \Y\ < 3'/ 5 the solid and dashed lines mark the edges and 
center of the slit, respectively. 
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Fig. 3. — (a) The kinematic integrals V as a function of the slit offset Y with respect to 
the major axis (Y = 0"). Open and filled circles correspond to slits crossing the secondary 
bar < 2'/8) and only the primary bar > 3'.'1), respectively, (b) The photometric 
integrals X as a function of the slit offset Y. (c) V as a function of X with different straight- 
line fits. These were obtained including the slits at Y — 0" and only the slits at \Y\ > 3"1 
(dashed line, slope f2 p sini = 8.0 ± 1.7 km s _1 arcsec^ 1 ) or the innermost slits at \Y\ < 2!8 
(solid line, slope tt Sjl sini = 63.7 ±7.1 km s _1 arcsec -1 ). The very different slopes of the two 
straight lines strongly suggest that the primary and secondary bars have different pattern 
speeds, (d) The residuals from the straight line fit to the slits at \Y\ > 3'.'1 and at Y — 0". 
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Fig. 4. — (a) The major-axis radial profile of the stellar line-of-sight velocity dispersion 
fitted with an exponential profile at r > 10" (solid line), (b) The major-axis radial profile 
of the stellar line-of-sight velocity (after subtracting the systemic velocity = 1312 ± 3 
km s _1 ) and the V c sini curve (solid line) with errors (dotted lines) obtained by applying the 
asymmetric drift for r > 30" as in ADC03. Dashed lines have slope Q St \ sini = 63.7 km s" 1 
arcsec -1 and 1 S;2 1 sin i = 15.1 km s" 1 arcsec" 1 . In (a) and (b) the measured profiles are 
folded around the center, with filled circles and asterisks referring to the SE (receding) and 
NW (approaching) sides, respectively. 



